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Abstract 

We stir with a focused laser beam a Bose-Einstein condensate of ®^Rb atoms confined 
in a magnetic trap. We observe the formation of a single vortex for a stirring frequency 
exceeding a critical value. At larger rotation frequencies we produce states of the condensate 
for which up to eleven vortices are simultaneously present. We present measurements of the 
decay of a vortex array once the stirring laser beam is removed. 

Pacs: 03.75.Fi, 67.40.Db, 32.80.Lg 

1 Introduction 

The discovery of Bose-Einstein condensation of atomic gases ||l], ^, |^, ^] has led to a new impulse 
in the physics of quantum gases. Among the several questions that can be studied in these 
systems, superfluidity is one of the most intriguing and fascinating. A first hint to the superfluid 
behaviour of these systems was provided by the study of the oscillation frequencies of the normal 
modes of these systems Q . A more direct evidence has been found by recent experiments aiming 
to study the energy deposited in the condensate by a moving "object" {i.e. the hole created by 
a blue detuned laser) @ and the rotational properties of the gas 0, ^. 

The relation between the rotational and superfluid properties of a gas or a liquid is illustrated 
by the famous "rotating bucket" experiment (see Fig. ||). When an ordinary fluid is placed in 
a rotating container, the steady state corresponds to a rotation of the fluid as a whole together 
with the vessel. The velocity field for a normal fluid is the same for a rigid rotating body: 

V = Q X r . (1) 



Superfluidity, first observed in liquid HeU, changes dramatically this behaviour ||10|, |12| [l 



14 1 . Indeed, for a low enough rotation frequency, the superfluid is not entrained at all by the 



walls of the rotating bucket, and it stays at rest in the laboratory frame. 

For a more quantitative description of the superfluid phenomenon, consider a Bose-Einstein 



condensate with weak, repulsive interactions (see e.g. [12, 13, |15[). This model is relevant for 
the ^'^Rb gas that we investigate experimentally since (i) the scattering length a characterizing 
the two-body interactions is positive, (ii) the "gaseous" parameter (pa^)^/^ is much smaller than 
1 {p is the average spatial density of the condensate). When cooled well below the condensation 
temperature, nearly all the atoms occupy the same state, and we can write the condensate wave 
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Figure 1: The rotating bucket experiment. A fluid is placed in a cylindrical container, which is rotated at a 
constant angular frequency CI along its axis. The velocity field at equilibrium reflects the superfluid character of 
the gas or the liquid. Indeed, for a low enough rotation frequency, the superfluid is not entrained at all by the 
walls of the rotating bucket, and it stays at rest in the laboratory frame. 



function as il^{r) = \/~p(r) e*'^*-^\ where p{r) is the spatial density. The velocity field is given 
by mv = h'VS, which is manifestly irrotational. This fact forbids, in particular, the rigid-body 
velocity field given in (|lj). An energy analysis, which is outlined in the next section of this paper 
for the particular case of harmonic confinement, leads to the following conclusions. For a low 
enough angular rotation frequency 17, the ground state is the same as that when the bucket 
is at rest. In this state, the phase of the condensate is constant over the whole volume and 
no rotational motion occurs. However above a critical frequency a line of singularity, i.e. 
a vortex line, appears. The condensate density drops to zero on axis, and the phase of the 
condensate wave function is now given by S{r) = 6, where 9 is the azimuthal angle around the 
vortex axis (oriented along z). The corresponding velocity field varies as r~^, where r is the 
distance to the vortex center. Moreover, the velocity field has a quantized circulation around 
the axis: 

2Trh 

vdr = , (2) 

m 

where m is the mass of the particles in the fluid. For frequencies notably higher than 17^, one 
might expect that a single vortex line with a quantum number n larger than 1 would appear 



{S{r) = n9). However this state is unstable |12, 13, 14| and fragments into n vortices each with 
a unit circulation quantum. 

These phenomenon have been observed in experiments with superfluid liquid helium; how- 
ever, the model of a dilute Bose gas presented in the previous paragraph is not applicable to 
this dense system. Let us quote here two milestones in this very rich field of research. In 1958 
Vinen performed an experiment in which he detected the presence of a single vortex filament 
by observing the induced frequency shift of the vibrational modes of a wire at the center of a 
rotating liquid helium bath [16|. In 1982 Yarmchuk and Packard obtained images of a vortex 
lattices in the superfluid by imaging electrons initially trapped at the cores of the vortex lines 



0. 
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Figure 2: A cigar-shaped atomic cloud confined in an axisymetric magnetic trap is stirred by a far-detuned laser 
beam. The laser beam propagates along the long axis of the cigar (z), and it creates a dipole potential which is 
anisotropic in the xy plane. This anisotropy rotates around the z axis at the angular frequency S7. 



2 Vortices in trapped gaseous condensates 



The generation of quantized vortices in trapped atomic gases has been the subject of numerous 
theoretical studies in the recent years. Two schemes have been considered. The first one uses 
laser beams to engineer the phase of the condensate wave function and produce the desired 
velocity field [|l8|, 19, 2C, 21, [2^. Recently this scheme ||2l| has been successfully applied to a 



binary mixture of condensates, resulting in a quantized rotation of one of the two components 
around the second one [^. Phase imprinting has also been used for the generation of solitons 
inside a condensate [23 




The second scheme, which is explored in the present work, is directly analogous to the 
rotating bucket experiment [25, |2^. The atoms are confined in a static, cylindrically-symmetric 
loffe-Pritchard magnetic trap upon which we superimpose a non-axisymmetric, attractive dipole 
potential created by a stirring laser beam (see figure ^). The combined potential leads to a cigar- 
shaped harmonic trap with a slightly anisotropic transverse profile. The transverse anisotropy is 
rotated at angular frequency 17 as the gas is evaporatively cooled to Bose-Einstein condensation, 
and it plays the role of the bucket wall roughness. 

In this scheme, the formation of vortices is a consequence of thermal equilibrium. In the 
frame rotating at the same frequency as the anisotropy, the Hamiltonian is time-independent and 
one can use a standard thermodynamics approach to determine the steady-state of the system. 
In this frame, the Hamiltonian can be written H = H — ilL^, where H is the Hamiltonian in the 
absence of rotation, and is the total orbital angular momentum along the rotation axis. For 
a gas with repulsive interactions, the term —QLz may favor the creation of a state where the 
condensate wave function has an angular momentum h along the z axis and therefore contains 
a vortex filament [ 



25, 29 



31 



33, 34, 35, 36 



This is illustrated in figure |^, which displays the lowest eigenenergies for the Gross-Pitaevski 
equation describing a gas with N atoms confined in a two-dimensional isotropic harmonic po- 
tential of frequency ut- 



Aw H — uiiot 

2m ^ 2 * 



ATrh'^Na , 



m 



(3) 



Here a is the scattering length characterizing the 2-body interaction. In figure |3|a, we consider 
an ideal gas (a = 0). The ground state found for = 0, which has zero angular momentum, 
remains the ground state of the system until the stirring frequency 0, reaches cot- Above this 
value, the gas is no longer stable since the centrifugal force mQ'^r exceeds the restoring force of 
the confining potential —mujfr. For a gas with repulsive interactions (figure |3|b), the energy of 
the state with angular momentum £ = 1 crosses the energy of the 1 = state at a critical value 
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Figure 3: Lowest transverse energy levels for the Gross-Pitaevski equation of a gas confined in a bi- dimensional 
harmonic potential with frequency uit, and stirred at a angular frequency n. (a) In absence of interaction the 
ground state of the system is the same as that for f2 = for any stirring frequency in the stability domain Q < uit- 
(b) For an interacting gas with repulsive interactions (assumed here to be weak), there exists a critical frequency 
Qc above which the ground state of the system has an angular momentum i = 1. The effect of the stirring 
anisotropy on the energy levels has been neglected. 

Qc < Indeed one can infer that for Q = u;± the energy of the state £ = 1 is strictly below the 
energy of the state i = 0: (i) these two states have the same energy in absence of interactions 
for 0, = uJt] (ii) repulsive interactions increase the energy of these states, by an amount which is 
larger for the 1 = state, since it has a smaller volume than the state with 1 = 1. 

When the interactions are large enough to reach the Thomas-Fermi limit, for which the 
kinetic energy is negligibly small compared to the potential and interaction energies, the solution 
of for i = 1 leads to a radius for the vortex core of the order of the healing length ^ = 
(Svrap)"^/^, where p is the density of the condensate at the center of the trap in the absence of 
a vortex . 

3 Experimental setup and determination of Qc 

Our set-up has been described in detail previously [^, ^, and we only briefly outline the main 
elements. The atoms are confined in an loffe-Pritchard magnetic trap. The slow oscillation 
frequency of the elongated magnetic trap is uJz/{2'7t) = 11.7 Hz (z is horizontal in our setup), 
while the transverse oscillation frequency is uj±/(27r) = 219 Hz. Bose-Einstein condensation is 
reached using evaporative cooling produced by a chirped radio-frequency source. For the data 
presented here, the final frequency of the evaporation ramp is chosen only a few kHz above the 
value which completely empties the trap. This garantees that the temperature of the atomic 
gas is below 80 nK, the uncondensed cloud being only marginally visible. The number of atoms 
in the condensate is of the order of 140 000. 

The stirring laser beam is switched on approximately at the time when the cloud reaches 
the condensation point. It propagates along the slow axis of the magnetic trap. The beam waist 
is Ws = 20.0 (lb 1) jum and the laser power P is 0.4 mW. The recoil heating induced by this 
far-detuned beam (wavelength 852 nm) is negligible. Two crossed acousto-optic modulators, 
combined with a proper imaging system, allow for an arbitrary translation of the laser beam 
axis with respect to the symmetry axis of the condensate. 

The motion of the stirring beam consists in the superposition of a fast and a slow compo- 
nent. The optical spoon's axis is toggled at a high frequency (100 kHz) between two symmetric 
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Figure 4: Absorption images of a Bose-Einstein condensate stirred with a laser beam (after a 27 ms time-of- 
fiight). For both images, the condensate number is A^o = (1-4 ± 0.5) 10^ and the temperature is below 80 nK. 
The rotation frequency f2/(27r) is respectively (a) 145 Hz; (b) 152 Hz. 



positions about the long trap axis z. The intersections of the stirring beam axis and the z = 
plane are ±a(cos0Ua;. + sin6'uj;), where the distance a is 8 //m. The fast toggle frequency is 
chosen to be much larger than the magnetic trap frequencies so that the atoms experience an 
effective two-beam, time-averaged potential. The slow component of the motion is a uniform 
rotation of the angle Q = OA,. The value of the angular frequency is maintained fixed during 
the evaporation at a value chosen between and 27r x 220 rad s~^. 

The dipole potential, proportional to the power of the stirring beam, is well approximated 
by vrnx)\(ex^'^ + ey^^)/2. The (X^Y) basis is rotated with respect to the fixed axes by 
the angle and ex = 0.03 and ey = 0.09 for the parameters given above [^. The action 
of this beam is essentially a slight modification of the transverse frequencies of the magnetic 
trap while the longitudinal frequency is nearly unchanged. The overall stability of the stirring 
beam on the condensate appears to be a crucial element for the success of the experiment, and 
we estimate that our stirring beam axis is fixed to and stable on the condensate axis to within 
2 fj,m. 

After the end of the evaporation ramp, we let the system reach thermal equilibrium in this 
"rotating bucket" for a duration t,- = 500 ms. The vortices induced in the condensate by the 
optical spoon are then studied using a time-of-flight analysis (27 ms) after the atoms have been 
released from the magnetic trap. Due to the atomic mean field energy, the initial cigar shape of 
the atomic cloud transforms into a pancake shape during the free fall. The transverse xy and 
longitudinal z sizes grow by a factor of 40 and 1.2 respectively p9| . In addition, the core size 
of the vortex should expand at least as fast as the transverse size of the condensate 4^, 41]. 



Therefore a vortex with an initial diameter 2^ = 0.4 ^m for our experimental parameters is 
expected to grow to a size of 16 ^m. 

At the end of the time-of-flight period, we illuminate the atomic sample with a resonant 
probe laser for 20 fis. The shadow of the atomic cloud in the probe beam is imaged onto a CCD 
camera with an optical resolution 7 fim. The probe laser propagates along the z-axis so that 
the image reveals the column density of the cloud after expansion along the stirring axis. 

For the trap parameters given above, we have found that no singularity appears in the cloud 
image as long as the stirring frequency is smaller than 149 Hz. However, for a stirring frequency 
equal to or above 150 Hz (but smaller than 160 Hz, see next section) a density dip systematically 
appears at the center of the cloud. Figure ^ illustrates this transition by showing two pictures 
of the condensate after expansion; the first with a stirring frequency O = 145 Hz (FigJ^) below 
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and the second Q = 152 Hz (Fig.p3) above the critical frequency. 



The dip in the optical density can reach up to 50 % of the maximal column density and 
it constitutes an unambiguous signature of the presence of a vortex filament. In particular it 
cannot be induced by a mechanical action of the stirring laser, which in fact creates a restoring 
force towards the center of the magnetic trap since it is detuned red of resonance. The fact that 
the density does not vanish in the dip may be due to several reasons: (i) The vortex filament 
may not be perfectly straight but rather an oscillating line (the so-called Thomson mode |p^ ). 
(ii) Some non-condensed atoms may be trapped in the core of the vortex, (iii) The core diameter 
after expansion of 18 fim at the half max of the dip is not much larger than the resolution 
limit of the imaging optics. 

Finally, we note that the critical frequency is notably larger than the predicted value of 



91 Hz (see also |27, 29, 34, 32|). This deviation may be due to the marginality of the 



Thomas- Fermi approximation for our relatively low condensate number. 



4 Dynamics of a vortex array 



When the stirring frequency is increased notably above the value i^c, multiple dips appear in 
the absorption image of the cloud. Each of these dips has approximately the same spatial width 
as the single dip of Fig. ^. For the trap parameters given above, the threshold frequency for 
the appearance for such multi-vortex patterns is 160 Hz, and condensates with up to 4 vortices 
have been observed |p. 

We have found that we can generate patterns with more than 4 vortices if we reduce the 
transverse trapping frequency. This can easily be done by adjusting the value of the bias 
magnetic field at the center of the trap. The data that we give now have been obtained with a 
transverse frequency of LOt / 2it = 169 Hz and a longitudinal frequency of tOz / 2Tr = 11.7 Hz, 
the same as before. 



Figure 5: Arrays of vortices in a Bose-Einstein condensate stirred by a laser beam (after a 27 ms time-of-flight) . 
These pictures have been obtained with a magnetic trap less stiff than for the pictures of Fig.^ {cot 2iy = 169 Hz). 

For an atom number similar to the previous data (~ 150 000), the critical frequency is now 
Q,c / 2tt = 119 Hz. This still corresponds to a ratio ^d^t ~ 0.7. When we increase the stirring 
frequency above 130 Hz, arrays of several vortices can be observed. For instance, Fig.^ shows 
images of condensates with 7, 8 and 11 vortices. These three images have all been obtained 
with the same experimental conditions [VL / 2tt = 135 Hz), apart from residual shot-to-shot 
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fluctuations. In the pictures with 7 and 8 vortices (Figj5|a and Fig.^), one vortex is found 
at the center of the trap, and the remaining ones form a regular hexagon or heptagon. The 
11-vortex pattern (Fig.||c) has 3 vortices forming an equilateral triangle surrounded by a regular 
octogon. Some of the dips of the octogon are hardly visible, since they are located at the border 
of the condensate cloud. 

In ^] we studied the lifetime of the one vortex pattern of Fig.^ when the stirring laser is 
removed. We found that this structure may survive for up to ~ 1 s before decaying, proving 
the metastability of the current associated with the vortex |42]. We now present similar data 



obtained for a multiple vortex structure. For this experiment we chose a rotation frequency of 
= 135 Hz and adjusted the number of atoms in the condensation to maximize the probability 
that a 5-vortex pattern form. For each point in the decay study, a vortex lattice is formed, then 
the stirring beam is turned off, and the condensate is allowed to evolve in the pure magnetic 
trap for an adjustable time. A time-of-flight analysis is performed to determine the number of 
vortices still present. 
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Figure 6: Decay of a vortex array, (a) Average number of vortices as a function of time spent by the gas in the 
axisymetric trap after the end of the stirring phase, (b) Fraction of the images displaying n vortices as a function 
of time. After 2 seconds all pictures exhibit a condensate with no vortex. 

The results are displayed in Fig.^ which presents the decay of the average number of vortices 
as a function of time, each point representing the average of 10 shots. The characteristic time 
for the reduction by a factor 2 of this average number is 750 ms. A more detailed representation 
is shown in Fig.p3 which displays the fraction of pictures showing n vortices (n = 0, . . . ,6) as 
a function of time. It seems apparent that the vortex pattern decays by losing one vortex at a 
time. For instance, at time t = 100 ms, 9 patterns out of 10 exhibit 5 vortices; 400 ms later, 
7 patterns out of 10 exhibit 4 vortices. It is remarkable that the vortex pattern readjusts itself 
shortly after the loss of a vortex. For instance, the 3-vortex patterns are found most often as 
a quasi-centered and quasi-equilateral triangle, while the 4-vortex patterns have essentially a 
square shape. 



5 Conclusion 

We have reported the formation of single and multiple-vortex structures in a gaseous Bose- 
Einstein condensate when it is stirred by a laser beam which produces a slight rotating 
anisotropy. We have also presented measurements of the lifetime of the pattern of vortices 
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when the rotating anisotropy is removed. 



Several extensions of this work can be considered. Direct evidence for the quantized circu- 
lation of the velocity field (Eq.Q) could be obtained either from an interferometric measurement 
similar to the one performed in @ (see e.g. |34] or ||2^), or from a study of the elementary 



excitations of the vortex filament p3|, 44, Hq]. Also the role of the thermal component in the 



nucleation and decay of the vortex pattern remains to be elucidated |32, 34, 47, 48 1. The 
rotational properties of the gas before the condensation point, where the collisional dynamics 
can be described by the classical or quantum Boltzmann equation, constitute other very inter- 
esting questions |4^. Finally, the scattering of quasi-resonant light by this medium may lead 
to spectacular phenomena, such as a "black hole" type behaviour when the speed of sound or 



mass flow in the condensate exceeds that of light |50|. 



We wish to dedicate this paper to the memory of Dan Walls who has had a seminal influence 
on the field of quantum optics in general, and in the understanding of the physics of gaseous 
Bose Einstein condensates in particular. 
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